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There is now abundant evidence for the presence of large quantities of unseen 
matter surrounding normal galaxies, including our own^'^. The nature of this 'dark 
matter' is unknown, except that it cannot be made of normal stars, dust, or gas, 
as they would be easily detected. Exotic particles such as axions, massive neutrinos 
or other weakly interacting massive particles (collectively known as WIMPs) have 
been proposed) but have yet to be detected. A less exotic alternative is normal 
matter in the form of bodies with masses ranging from that of a large planet to a few 
M0. Such objects, known collectively as massive compact halo objects^ (MACHOs) 
might be brown dwarfs or 'Jupiters' (bodies too small to produce their own energy 
by fusion), neutron stars, old white dwarfs, or black holes. Paczynski^ suggested that 
MACHOs might act as gravitational microlenses, occasionally causing the apparent 
brightness of distant background stars temporarily to increase. We are conducting a 
microlensing experiment to determine whether the dark matter halo of our galaxy is 
made up of MACHOs. Here we report a candidate for a microlensing event, detected 
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by monitoring the light curves of 1.8 milUon stars in the Large Magellanic Cloud for 
one year. The light curve shows no variation for most of the yesir of data taking, and 
an upward excursion lasting over 1 month, with a maximum increase of 2 mag. The 
most probable lens mass, inferred from the duration of the event, is 0.1 Mq. 

The MACHO Project'''* uses the gravitational microlcns signature to search for evidence of 
MACHOs in the Galactic halo, which is thought to be at least three times more massive than 
the visible disk^. (Two other groups are attempting a similar search^' ^°.) If most of our Galaxy's 
dark matter resides in MACHOs, the 'optical depth' for microlensing towards the Large Magellanic 
Cloud (LMC) is about 5 x 10~^ (independent of the mass function of MACHOs) , so that at any given 
time about one star in two million will be microlcnsed with an amplification factor A > 1.34 (rcf. 5). 
Our survey takes advantage of the transverse motion of MACHOs relative to the line-of-sight from 
the observer to a background star. This motion causes a transient, time-symmetric and achromatic 
brightening that is quite unlike any known variable star phenomena, with a characteristic timescale 
t = 2rE/v± where Ve is the Einstein ring radius and v± is the MACHO vdocil x- transverse to the 
linc-of-sight. For typical halo models, the time with A > 1.34 lOO-^Mmacho/M© days^ (where 
Mq is the mass of the sun). The amplification can be large, but these events are extremely rare; 
for this reason our survey was designed to follow > ten million stars over several years. 

The survey employs a dedicated 1.27m telescope at Mount Stromlo. A field-of-view of 0.5 
square degrees is achieved by operating at the prime focus. The optics include a dichroic beam- 
splitter which allows simultaneous imaging in a 'red' beam (6300 — 7600A) and a 'blue' beam 
(4500 — 6300A). Two large charge coupled device (CCD) cameras^^ are employed at the two foci; 
each contain a 2 x 2 mosaic of 2048 x 2048 pixel Loral CCD imagers. The 15 /im pixel size cor- 
responds to 0.63 arcsec on the sky. The images are read out through a 16 channel system, and 
written into dual ported memory in the data acquisition computer. Our primary target stars are 
in the LMC. We also monitor stars in the Galactic bulge and the Small Magellanic Cloud. As of 
15 September 1993, over 12000 images have been taken with the system. 

The data are reduced with a crowded-field photometry routine known as Sodophot, derived 
from Dophot^^. First, one image of each field that was obtained in good seeing is reduced in 
a manner similar to Dophot to produce a 'template' catalog of star positions and magnitudes. 
Normally, bright stars are matched with the template and used to determine an analytic point 
spread function (PSF) and a coordinate transformation. Photometric fitting is then performed on 
each template star in descending order of brightness, with the PSF for all other stars subtracted 
from the frame. When a star is found to vary significantly, it and its neighbors undergo a second 
iteration of fitting. The output consists of magnitudes and errors for the two colors, and six 
additional useful parameters (such as the of the PSF fit and crowding information). These are 
used to fiag questionable measurements that arise from cosmic ray events in the CCDs, bad pixels, 
and so on. 

These photometric data are subjected to an automatic time-series analysis which uses a set 
of optimal filters to search for microlensing candidates and variable stars (which we have detected 
in abundance^^). For each microlensing candidate a light-curve is fitted, and the final selection is 
done automatically using criteria (for example signal-to-noise, goodness of fit, achromaticity, color) 
that were established empirically using Monte Carlo addition of fake events into real light-curves. 

This analysis has been done on four fields near the center of the LMC, containing 1.8 million 
stars, with approximately 250 observations for each star. The candidate event reported here occurs 
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in the light-curve of a star at coordinates a = 05 14 44.5,(5 = —68 48 00 (J. 2000). (A finding 
chart is available on request from C. A.). The star has median magnitudes V ~ 19.6, R ~ 19.0, 
consistent with a clump giant (metal-rich helium core burning star) in the LMC. These magnitudes 
arc estimated using color transformations from our filters to V and R that have been derived from 
observations of standard stars. 

Our photometry for this star, from July 1992 to July 1993 is shown in Figure 1, and the 
candidate event is shown on an expanded scale in Figure 2, along with the color light curve. The 
color changes by < 0.1 mag as the star brightens and fades. A mosaic, showing portions of some 
of the CCD images used, is shown in Figure 3, with the relevant star at the center. The integrated 
number of PSF photoelectrons detected above the sky background in the template image is ~ 10*^, 
for a 300 sec exposure. The increase in counts during the peak is highly significant, as is clear from 
the Figures. 

Also shown in Figure 2 is a fit to the theoretical microlensing light-curve (see ref. 6) . The four 
parameters fit are (1) the baseline flux; (2) the maximum amplification A^ax = 6.86 ± 0.11; (3) 
the duration i = 33.9 ± 0.26 days; (4) the centroid in time 433.55 ± 0.04 days. The quoted errors 
are formal fit errors. Using the PSF fit uncertainties as determined by the photometry program, 
the best-fit microlensing curve gives a pcr degree of freedom of 1.6 (for 443 d.o.f.). 

A number of features of the candidate event are consistent with gravitational microlensing: the 
light curve is achromatic within measurement error, and it has the expected symmetrical shape. If 
this is a genuine microlensing event, the mass of the deflector can be estimated. Since the duration 
depends upon the lens mass, the relative velocity transverse to the line-of-sight, and the distance 
to the lens, (none of which are known), the lens mass cannot be uniquely determined from the 
duration. However, using a model of the mass and velocity distributions of halo dark matter, one 
can find the relative probability that a MACHO of mass m gave rise to the event. Thus, if this is 
geniune microlensing. Fig. 9 of Ref. 5 implies the most likely mass is approximately O.I2M0, with 
masses of O.O3M0 and 0.5Mq being roughly half as likely. However, this method does not properly 
take into account our detection efficiencies, and should be considered only a rough estimate. 

The mass range given above includes brown dwarfs and main sequence stars. Any microlensing 
star is very unlikely to be a red dwarf of the Galactic stellar halo, because one can show that the 
optical depth r* for microlensing by main sequence stars of the stellar halo is very low. Even if 
the mass function of the stellar halo rises as steeply as dN/dM oc M~^, as suggested recently^^ 
(here N is the number of stars per unit stellar mass interval) , r* is still a few hundred times smaller 
than the 5 x 10~^ optical depth estimated for MACHO microlensing. The chance of finding such 
a stellar microlensing event among our 1.8 million stars is therefore very small. 

The prospects for direct observation of a lensing object are not favorable. Even a star of 
O.5M0, for example, would have V ~ 24, and for many years would be within a small fraction of 
an arcsecond of the much brighter LMC star. 

We emphasize that the observed stellar brightening could be due to some previously unknown 
source of intrinsic stellar variability. The fit discrepancy near the peak is not yet understood; a 
more refined analysis of the data is underway. We do not yet have a spectrum of the star. A 
crucial test of the hypothesis that we are seeing gravitational microlensing by MACHOs in the 
galactic halo will be the detection of other candidates. So far, we have analysed only ~ 15% of our 
first year's frames, and we plan to continue observations until 1996; this will allow us to determine 
whether or not gravitational microlensing is really the cause. Additional events should show the 
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theoretically expected distribution of maxima, and should be representative of both the color- 
magnitude diagram and the spatial structure of the LMC. No repeats should be seen in any given 
star! (While this paper was in preparation, we were informed by J. Rich (personal communication) 
of the candidate events reported by the EROS collaboration. Note that the two groups use different 
definitions of characteristic time.) 

If such candidates do result from microlensing we should be able to determine the contribution 
of MACHOs to the dark matter in the Galactic halo. The results presented here encourage us to 
believe this will happen. 

We are very grateful for the skilled support given our project by the technical staff at the 
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is supported in part by the Office of Science and Technology Centers of NSF under cooperative 
agreement AST-8809616. Work performed at MSSSO is supported by the Bilateral Science and 
Technology Program of the Australian Department of Industry, Technology and Commerce. KG 
acknowledges a DOE OJI grant, and CWS thanks the Sloan Foundation for their support. 
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Figure Captions 



Figure 1 



The observed light-curve with estimated ±lcr errors. The upper panel shows 
Af,iue ) the flux (in linear units) divided by the median observed flux, in the blue 
passband. The lower panel is the same for the red passband. The smooth curve 
shows the best-fit theoretical microlensing model, fitted simultaneously to both 
colors. Time is in days from 1992 Jan 2 00:00 UT. 



Figure 2 



As in Figure 1, with expanded scale around the event. The curve shows the 
theoretical microlensing model, fitted simultaneously to both colors. The bottom 
panel is the color light curve, showing the ratio of red to blue flux, normalized 
so that the median is unity. 



Figure 3 



Selected Red CCD frames centered on the microlens candidate, showing obser- 
vations before, during and after the event. 
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